Constraining the high-density behavior of nuclear equation of state from strangeness 

production in heavy-ion collisions 
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The dynamics of pions and strange particles in heavy-ion collisions in the region of 1A GeV 
energies is investigated by the lanzhou quantum molecular dynamics model for probing the nuclear 
equation of state at supra-saturation densities. The total multiplicities and the ratios obtained in 

197 Au+ 197 Au 

over 12 C+ 12 C systems are calculated for selected Skyrme parameters SkP, SLy6, Ska 
and SIII, which correspond to different modulus of incompressibility of symmetric nuclear matter and 
different cases of the stiffness of symmetry energy. A decrease trend of the excitation functions of the 
ratios for strange particle production with increasing incident energy was observed. The available 
data of K + production measured by KaoS collaboration are described well with the parameter SkP, 
which results in a soft equation of state. The conclusions can not be modified by an in-medium 
kaon-nucleon potential. 
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Particle production in relativistic heavy-ion collisions 
has been investigated as a useful tool to extract the infor- 
mation of the nuclear equation of state (EoS) under ex- 
treme conditions in terrestrial laboratories, such as high 
density, high temperature and large isospin asymmetry 
etc. Besides nucleonic observables such as rapidity dis- 
tribution and flow of free nucleons or light clusters (such 
as deuteron, triton and alpha etc.), also mesons emitted 
from the reaction zone can be probes of the hot and dense 
nuclear matter. In particular, the strangeness produc- 
tions (K°' + , A and £~' 0,+ ) in heavy-ion collisions in the 
region of 1A GeV energies have been extensively investi- 
gated both experimentally [U-EI] and theoretically 0-0] ■ 
Kaons (K° and K + ) as a probe of EoS are produced in 
the high density phase without subsequent reabsorption 
process. It was noticed that the kaon yields are sensi- 
tive to the EoS in theoretical investigations by transport 
models. The available data already favored a soft EoS 
at high densities. The K°/K + ratio was also proposed 
as a sensitive probe to constrain the high-density behav- 
ior of the symmetry energy 0, . It is not only in un- 
derstanding the reaction dynamics, the high-density be- 
havior of the EoS also has an important application in 
astrophysics, such as the structure of neutron star, the 
cooling of protoneutron stars, the nucleosynthesis during 
supernova explosion of massive stars etc jjil EH ■ 

The lanzhou quantum molecular dynamics (LQMD) 
model has been successfully applied to treat the dynamics 
in heavy-ion fusion reactions near Coulomb barrier and 
also to describe the capture of two heavy colliding nu- 



clides to form a superheavy nucleus [12i . 1 1 31 ] . Further im- 
provements of the LQMD model have been performed in 
order to investigate the dynamics of pion and strangeness 
productions in heavy-ion collisions and to extract the in- 
formation of isospin asymmetric EoS at supra-saturation 
densities 14 - Tl7j j . In this work, the high-density behavior 
of nuclear equation of state is to be investigated through 
particle productions in heavy-ion collisions in the region 



of 1A GeV energies. 

In the LQMD model, the time evolutions of the 
baryons and mesons in the system under the self- 
consistently generated mean-field are governed by Hamil- 
ton's equations of motion, which read as 

OH . _dH 

Pi n j r i — T. • 

OTi dpi 



(1) 



Here we omit the shell correction part in the Hamiltonian 
H as described in Ref . [3] • The Hamiltonian of baryons 
consists of the relativistic energy, the effective interaction 
potential energy and the momentum dependent part as 
follows: 



H B = ^ V Pi + m i + Uint + U ™°r> 



(2) 



Here the p i and mi represent the momentum and the 
mass of the baryons. The momentum dependent term is 
taken as the same form in Ref. jl8j, which reduces the 
effective mass in nuclear medium. 

The effective interaction potential is composed of the 
Coulomb interaction and the local interaction 

Uint = Ucoul + Ui oc . (3) 

The Coulomb interaction potential is calculated by 

(4) 



Ucoul ^er/(r«/>/4L) 



where e is the charged number including protons and 
charged resonances. The = |rj — r j | is the relative 
distance of two charged particles. 

The local interaction potential energy is derived di- 
rectly from the Skyrme energy-density functional and ex- 
pressed as 



(5) 



Uioc = / Vi oc {p(r))dr. 
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The local potential energy-density functional reads 
„2 a 1+7 



V loc (p) 



2 Po 1 + 7 2/30 



o sur 
_ 2p 



[V(p n - ft,)] 5 



P" P" ' ' P 8/3 \ ,2 



+9rP S/3 / P l /3 , 



(6) 



where the p n , p p and p = p„ + p p are the neutron, 
proton and total densities, respectively, and the S = 
(pn ~ Pp)/{Pn + Pp) is the isospin asymmetry. Here, all 
the terms in the Skyrme energy functional are included 
in the model besides the spin-orbit coupling. The coef- 
ficients a, /3, 7, g sur , g\ S u r , 9t are related to the Skyrme 
parameters to,ti,t2,ts and Xq,X\,X2,x^ as [l3j |. 
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^ = ^(9*1-5^-4.^2)^, 



(7) 
(8) 



= ~il^( 2x i + + ^(2^2 + l)]po, (9) 
2 64 

5r = ^ (|tt 2 ) 7 (3t x + 5t 2 + Ax 2 t 2 )pf\ (10) 

The parameters of the potential part in the bulk sym- 
metry energy term are also derived directly from Skyrme 
energy-density parameters as 

a sy m = -g(2xo + l)*oPo, b sym = — — (2x 3 + l)*3Po> 
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= (J^J Po /3 [3*i^i - t 2 (5.x 2 + 4)]. (11) 

In the calculations, we use the Skyrme parameters SkP, 
SLy6, Ska and SIII, which result in different modulus of 
the incompressibility of symmetric nuclear matter and 
also different stiffness of the symmetry energy. In Table 
1 we list the LQMD parameters related to several typical 
Skyrme forces. 

Analogously to baryons, the evolution of mesons (here 
mainly pions and kaons) is also determined by the Hamil- 
tonian, which is given by 

N M 

HM='£(vf oul +u>(p i ,p i )). ( 12 ) 

i=l 

Here the Coulomb interaction is given by 

N B 

K Coul = V^, (13) 



E, 

3 = 1 y 
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where the Am and Nb are the total numbers of mesons 
and baryons including charged resonances. Here, we 



use the energy in vacuum for pions, namely tvlp^pi) = 
ypf + m 2 M with the momentum p t and the mass rau of 
the pions. We consider two scenarios for kaon (antikaon) 
propagation in nuclear medium, one with and one with- 
out medium modification. From the chiral Lagrangian 
the kaon and antikaon energy in the nuclear medium can 
be written as 0, [lj| 

W (P 4 , Pi) = [ m K + Pi - a Kpf + {bx Pi) 2 } 1/2 + b K Pt, 

(14) 



u(Pi,Pi) 



K 



+ P 2 - a-wpf + {bxPif 
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bxPt, 



(15) 

respectively. Here the b K = 3/(8/ 2 ) w0.32 GeVfm 3 , 
the clk and a-^ are 0.18 GeV 2 fm 3 and 0.3 GeV 2 fm 3 , re- 
spectively, which result in the strength of repulsive kaon- 
nucleon potential and of attractive antikaon-nucleon po- 
tential with the values of 25.5 MeV and -96.8 MeV at 
saturation baryon density. 

A hard core scattering in two particle collisions is as- 
sumed in the simulation of the collision processes by 
Monte Carlo procedures, in which the scattering of two 
particles is determined by a geometrical minimum dis- 
tance criterion d < yjQ.latot Jjf fm weighted by the Pauli 
blocking of the final states [2fJ, |21 1 . Here, the total cross 
section a t ot in mb is the sum of the elastic and all inelas- 
tic cross section. The probability reaching a channel in a 
collision is calculated by its contribution of the channel 
cross section to the total cross section as P c h = a c h / <J to t ■ 
The choice of the channel is done randomly by the weight 
of the probability. The primary products in nucleon- 
nucleon (NN) collisions at the 1A GeV energies are the 
resonances A(1232), N*(1440) and the pions. The reac- 
tion channels are given as follows: 

AA^AA, jvjVoM*, AAoAA, 

A o Att, A* O Att, AA -+ AAtt(s - state). (16) 

The cross sections of each channel to produce resonances 
are parameterized by fitting the data calculated with the 
one-boson exchange model [22}. In the 1 A GeV region, 
there are mostly A resonances which disintegrate into a 
7r and a nucleon, however, the A* yet gives considerable 
contribution to the high energetic pion yield. The en- 
ergy and momentum dependent decay width is used in 
the calculation [TjJ ■ The strangeness is created by inelas- 
tic hadron-hadron collisions. We include the channels as 
follows: 

BB -+ BYK, BB -+ BBKK, Bit -> YK, Bit -> NKK, 
Ytt -> AA, AA -> Yn, YN -> AAA. (17) 

Here the B strands for (N, A) and Y(A, S), K(K°, K+) 
and K(K°, K~). The parameterized cross sections of 
each isospin channels BB — » BYK [23| are used in the 
calculation. We take the parametrizations of the chan- 
nels Bit — > YK [24j besides the Ntt —> AK reaction 
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TABLE I. LQMD parameters and properties of symmetric nuclear matter for Skyrme effective interactions after the inclusion 
of the momentum dependent interaction. 



Parameters 


SkM* 


Ska 


SIII 


SVI 


SkP 


RATP 


SLy6 


a (MeV) 


-325.1 


-179.3 


-128.1 


-123.0 


-357.7 


-250.3 


-296.7 


/3 (MeV) 


238.3 


71.9 


42.2 


51.6 


286.3 


149.6 


199.3 


7 


1.14 


1.35 


2.14 


2.14 


1.15 


1.19 


1.14 


s „ r (MeV fm 2 ) 


21.8 


26.5 


18.3 


14.1 


19.5 


25.6 


22.9 


r r (MeV fm 2 ) 


-5.5 


-7.9 


-4.9 


-3.0 


-11.3 


0.0 


-2.7 


9r (MeV) 


5.9 


13.9 


6.4 


1.1 


0.0 


11.0 


9.9 


E sym (MeV) 


30.1 


33.0 


28.2 


27.0 


30.9 


29.3 


32.0 


a sym (MeV) 


62.4 


29.8 


38.9 


42.9 


94.0 


79.3 


130.6 


b sym (MeV) 


-38.3 


-5.9 


-18.4 


-22.0 


-63.5 


-58.2 


-123.7 


Csym (MeV) 


-6.4 


-3.0 


-3.8 


-5.5 


-13.0 


-4.1 


12.8 


poo (fm -3 ) 


0.16 


0.155 


0.145 


0.144 


0.162 


0.16 


0.16 


Koo (MeV) 


215 


262 


353 


366 


200 


239 


230 



[25j | . The results are close to the experimental data at 
near threshold energies. The cross section of antikaon 
production in inelastic hadron-hadron collisions is taken 
as the same form of the parametrization used in hadron 
string dynamics (HSD) calculations J26]. Furthermore, 
the elastic channels are considered through the channels 
KB — > KB and KB —t KB and we use the parametriza- 
tions in Ref. [271 ]. 

The dynamics of pion and strangeness production in 
heavy-ion collisions in the region of 1A GeV energies is 
investigated systematically by using the LQMD model. 
Shown in Fig. 1 is a comparison of charged pions and 
strange particles in the 197 Au+ 197 Au reactions and com- 
pared with the available data measured by FOPI col- 
laboration 28]. The charged pions are overpredicted at 
near threshold energies, in particular for the 7r - produc- 
tion, which results from the different stiffness of the sym- 
metry energy [TBI 16 1. A similar trend is clear in the 
strangeness production at different Skyrme forces. Over- 
all, the parameter SIII leads to lower yields of strange 
particles. In order to reduce some uncertainties, such as 
cross sections, in-medium decay width of resonances etc., 
we calculated the double ratios of charged pions per mass 
numbers of reaction partners produced in 197 Au+ 197 Au 
over 12 C+ 12 C collisions as shown in Fig. 2. A slight 
increase of the value with incident energy is obvious for 
each parameters, which is consistent with the experimen- 
tal observations 0] ■ One can see that the influence of the 
equation of state on the excitation functions is slightly. 
Therefore, it is difficulty to extract the high-density infor- 
mation of nuclear equation of state from pion emissions. 

Neutral strange particles can not be influenced by the 
Coulomb interaction of surrounding baryons. Production 
of A is of great interest as a probe of high-density infor- 
mation of nuclear matter. We calculated the excitation 
functions of double ratios of A emissions in heavy-ion 
collisions as shown in Fig. 3. Different distributions are 
pronounced for each Skyrme forces. The parameters Sly6 
and SkP give a monotonically decreasing trend with in- 



creasing incident energy. However, the forces SIII and 
Ska reduce the values and even slightly drop at low en- 
ergy, which correspond to a hard equation of state of 
symmetric nuclear matter. Experimental measurements 
are helpful for constraining the stiffness of nuclear EoS. 
The increase of the ratios with decreasing incident en- 
ergy is caused by the fact that the number of collisions 
where the involved particles are encountered prior to the 
production of a kaon does not reach a sharp region for 
the system C+C at far low energies, which is related to 
higher compressible nuclear matter. 

To extract more information of the high-density EoS, 
in Fig. 4 we show the kaon production in the 197 Au+ 
197 Au and 12 C+ 12 C reactions for head-on collisions, nor- 
malized by the corresponding mass numbers and com- 
pared with the KaoS data [2|. A similar distribution 
with the A emission in Fig. 3 is observed, which re- 
sults from the same channels for the kaon and A pro- 
duction, such as BB — > BAK and Bit — > AK. A larger 
high-density region (> po) is formed in 197 Au+ 197 Au 
collisions and the compression depends on the nuclear 
equation of state. Whereas the compression in 12 C+ 12 C 
reaction is small and not sensitive to the stiffness of the 
EoS. The comparison to the available data further sup- 
ports a soft equation of state in the high-density region, 
in particular the parameter SkP, where the increase of 
the double ratios with decreasing incident energy is more 
pronounced when going far below threshold energies. In- 
clusion of the in-medium kaon-nucleon (KN) potential as 
shown in Fig. 5 almost does not change the results. The 
conclusions are consistent with the calculations by Fuchs 
et al. within the QMD model [6] and by Hartnack et 
al. within the Isospin-QMD (IQMD) model 0- Figure 6 
is the incompressibility dependence of the double ratios 
at incident energies 0.8A GeV and 1.2A GeV. The KN 
potential slightly reduces the ratios because of the de- 
crease of the secondary collision rate for the production 
of kaon. Furthermore, the soft equation of state is close 
to the experimental data. 
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FIG. 1. (Color online) Excitation functions of charged pions and strangeness production in central 197 Au+ 197 Au collisions at 
different Skyrme forces. 




FIG. 2. (Color online) Comparison of the calculated excitation functions of the ratios of charged pion multiplicities obtained 
in 

197 Au+ 197 Au 

over 12 C+ 12 C reactions for head-on collisions at different Skyrme forces. 



In summary, the production of pions and strange par- 
ticles in heavy-ion collisions in the region of 1A GeV 
energies as probes of the nuclear EoS is investigated 
within the framework of the LQMD model. The charged 
pions of the available experimental data can be repro- 
duced well within the selected Skyrme parameters SkP, 
SLy6, Ska and SIII. The excitation functions of the dou- 
ble ratios of charged pions obtained in the 197 Au+ 197 Au 
over 12 C+ 12 C collisions weakly depend on the nuclear 
equation of state. Whereas, the ratios of strange parti- 
cle production are sensitive to the potential parameters 
of nucleon-nucleon interaction. The comparison to the 
KaoS data favors a soft equation of state. This conclu- 
sion can not be modified by the input parameters, such 
as the KN potential, the lifetime of the A in the nuclear 



matter, cross sections in production of strangeness. 
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FIG. 4. (Color online) Excitation functions of kaon multiplicities calculated by the LQMD model without KN potential in 
central 197 Au+ 197 Au over 12 C+ 12 C collisions and compared with the available data from KaoS collaboration for K + production. 
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FIG. 5. (Color online) The same as in Fig. 4, but for calculations with KN potential. 



7 



2 



data range 



E, =0.8 A GeV 

lab 

o w/o KN pot. 
• with KN pot. 



200 



250 



300 



350 



data range 



s 1 



E, =1.2 A GeV 

lab 

o w/o KN pot. 
• with KN pot. 



200 



250 



300 



350 



FIG. 6. (Color online) Double ratios calculated within the 
LQMD model at 0.8A GeV (up panel) and 1.2A GeV (down 
panel) as a function of the incompressibility coefficient Ksm 
and compared with the available data for the K + production. 



